1

Import the ESRI shapefile of German districts and the district attribute table. Join the two data frames, transform the CRS to EPSG:3035 and check your changes.

You need to rename one of the id variables or adjust your join accordingly (AGS = district_id).

2

We want a first descriptive visual of the distribution of Covid-19 deaths in Mannheim and the surrounding districts. Calculate the number of Covid-19 deaths (death_rate) by population (population) and multiply with 100,000.

Filter the district of Mannheim (district_id == "08222"), find the surrounding districts, and plot Mannheim and its surrounding districts.

You can use the dplyr function sf::bind_rows() to combine the two spatial objects, “Mannheim” and “Mannheim Surroundings”.

3

Save your data set of Mannheim and its surrounding districts as an ESRI Shapefile.